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PHASE  SHIFT  ERRORS  ASSOCIATED  WITH  THE 
WKB  APPROXIMATION 

I.  INTRODUCTION 

In  tms  paper  a  numerical  technique  which  yields  an  exact  solution  to 
the  one-dimensional  scattering  problem  is  presented.  The  method  is  easy  to 
implement  in  computer  code  and  it  is  completely  general  in  that  no  special 
assumptions  concerning  the  scattering  potential  have  to  be  made.  Although 
tne  procedure  is  describee  in  terms  of  a  scattering  potential,  it  will  be 
shown  that  problems  dealing  with  electromagnetic  propagation  in  inhomogeneous, 
stratified  piasmas  or  dielectrics  which  are  described  in  terms  of  an  index  of 
refraction  can  be  reformulated  in  terms  of  potentials.  For  the  case  of  wave 
propagation  in  piasmas,  it  turns  out  the  associated  potential  is  proportional  • 
to  the  electron  density.  The  usual  method  which  is  applied  to  general 
scattering  problems  is  the  well  known  WKB  method  and  the  associated  connection 
formulas.  A  statement  for  the  validity  of  the  WKB  approach  is  that  thr*  change 
of  tne  local  wave  number  over  a  distance  equal  to  the  local  wavelength  must  be 
small  in  comparison  to  the  wave  number.  To  apply  the  connection  formulas  it 
is  assumed  that  in  the  neighborhood  of  a  point  where  the  local  wave  number  is 
zero  the  potential  can  be  represented  by  a  linear  function  with  a  non-zero 
slope.  There  are  many  interesting  problems  where  these  conditions  are  not 
satisfied.  Given  an  inhomogeneous  plasma  or  dielectric,  as  one  moves  to  lower 
frequency  electromagnetic  waves,  the  wavelengths  become  comparable  to  the 

scale  of  the  variation  of  media  properties.  If  the  frequency  of  the  wave  is 
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such  that  the  index  of  refraction  is  approximately  zero  over  a  significant 
path  length,  tnen  the  assumed  conditions  for  the  WKB  approximation  and  the 
connection  formulas  will  not  be  satisfied. 

The  orientation  of  this  paper  is  to  investigate  the  phase  error  associated 
with  the  WKB  approximation.  It  is  common  in  the  literature  to  introduce  an 
additional  phase  shift  of  TT/ 2  for  reflected  waves  over  and  above  the  phase 
shift  resulting  from  tne  geometric  optics  approximation  (see  [1]).  It  has 
been  found  that  this  approach  can  be  fairly  accurate  even  if  some  of  the 
assumed  conditions  of  Che  WKB  approximation  are  only  approximately  satisfied. 
Conversely,  there  are  also  interesting  problems  for  which  the  WKB  approximation 
introduces  a  significant  error  even  when  the  potential  is  slowly  varying 
relative  to  the  wavelength.  One  situation  for  which  tne  WKB  method  results  in 
large  phase  errors  is  a  "potential  IedgeM .  Ledge  type  potentials  can  often  be 
associated  with  HF  propagation  in  the  ionosphere,  since  the  daytime  index  of 
refraction  for  the  top  of  the  E-layer  is  frequently  approximately  constant 
for  distances  up  to  50  kilometers.  It  has  also  been  found  that  when  the 
scattering  potential  is  in  the  shape  of  a  well  for  which  there  are  resonance 
(metastable)  states  as  discussed  by  Bohm  [2],  then  the  phase  error  may  be 
quite  small  except  within  a  very  narrow  range  about  the  resonance  frequency 
and  within  this  narrow  band  the  phase  error  fluctates  rapidly.  For  the  case 
of  the  ionosphere,  a  potential  well  corresponds  to  a  local  minimum  in  the 
electron  density. 

The  problem  to  be  solved  is  that  of  a  plane  wave  incident  on  a  scattering 
potential.  For  the  corresponding  one  dimensional  problem,  the  potential  is 
non-zero  only  in  a  neighborhood  of  the  origin.  The  amplitudes  and  phases  of 
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tne  reflected  ana  transmitted  waves  are  described  by  the  complex  reflection 


coefncient  R  and  transmission  coefficient  T.  The  basic  situation  is  a 
two-sided  problem  with  energy  being  reflected  and  transmitted.  A  one-sided 
problem  can  be  defined  as  one  in  which  tne  scattering  potential  is  such  that 
tne  transmission  coefficient  is  zero.  The  technique  initially  developed  is 
restricted  to  two-sided  proolems  in  that  it  requires  the  existence  of  a  non¬ 
zero  transmission  coefficient.  however  a  one-sided  problem  can  be  converted 
into  a  two-siaed  prooiem  by  "truncating”  the  scattering  potential  on  the  side 
opposite  to  tne  incident  wave  in  such  a  way  that  the  transmission  coefficient 
is  very  small  ,  i.e.  j  T  j  ^  ~  lO"*1^.  In  fact,  numerical  results  have  shown 
tor  this  case  that  R  is  essentially  unchanged  as  long  as  the  truncation  point 
is  cnosen  such  that  j T |  ^ 

Since  the  scattering  of  a  wave  by  a  known  potential  is  described  by  a 
differential  equation  it  seems  natural  to  apply  a  numerical  scheme  such  as  the 
Runge-Kutta  method  to  obtain  an  exact  solution.  A  problem  arises  because  the 
reflection  and  transmission  coefficients  are  initially  unknown  and  thus  the 
initial  conditions  necessary  for  a  numerical  integration  are  not  available. 

It  is  shown  tnat  for  the  original  scattering  problem  there  is  an  associated 
problem  for  which  the  initial  conditions  on  the  far  side  of  the  scattering 
potential  can  be  specified.  By  solving  the  associated  problem  and  then 
transforming  back  to  the  original  problem,  one  obtains  R  and  T. 

II.  ELECTROMAGNETIC  PROPAGATION  'FORMULATED  IN  TERMS  OF  A  SCATTERING  POTENTIAL 

Assume  that  in  a  plane  stratified  medium,  the  index  of  refraction  is  given 
by  n<,x;  for  -  00  <  x  <  00  #  For  an  lsotropic  plasma,  the  expression  for  n(x) 
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is 


n^xj  =  i  -  XI. x) 

wnere  Xlx;  _>0.  For  a  dielectric  n(.x)  is  taken  to  be 

n2(.x;  =  i  +u2vxj 

with  ^.U.  Jioth  cases  can  be  described  by  writing 

n^(.x;  =  1  +  r (xj .  (1) 

Tne  wave  equation  tor  propagation  in  tnis  medium  has  the  form 

V2W(.x,y,z)  +  k^n-^(x)  W(x,y,zJ  =0.  _  (2) 

here  x  is  tne  wave  number  and  is  related  to  the  free  space  wavelength  (A.)  by 
x  =  2  it/  X  . 


The  incident  plane  wave  is  given  by 


WAx,y,z)  =  u (.x jexpl  ik Cycos  0  +  zcos  y  )]  (3) 

wnere  tne  direction  of  propagation  is  defined  by  the  direction  cosines  (cos  a, 
cos  g,  cos  7  ).  If  i.3)  and  (U  are  used  in  (2)  we  find 


(4) 
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wn  ere 


E  -  iC  2-  COS  2  CL 

V^x)  =  -k2r(x). 

Except  for  the  constraint  on  the  potential,  equation  (4)  is  the 
one-dimens  lonai  Scnrodmger  equation.  For  a  plasma  tne  scattering  potential 
satis  ties 


VCxJ  >0  all  x, 


wnile  for  a  dielectric 


V(x)  <  0  all  x. 


III.  SOLUTION  OF  THE  SCHROD LINGER  EQUATION 


The  equation  to  be  solved  is 

,2 


7? +  [ 


k  -V(.x)  .  u  =  0  • 


(5) 


For  simplicity  normal  incidence  of  the  plane  wave  is  assumed.  It  is  also 
assumed  tnat  Vlx;  >  0  and  VCx)  =  0  outside  the  interval  (0,L).  The  system  to 
be  solved  then  Decomes 
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u(x)  =  e1^*  +  x<  0 

(6) 

u(x;  =  Te^-R*  ,  x  >L 

/ 

where  R  and  T  are  the  unknown  reflection  and  transmission  coefficients.  A 
system  which  is  equivalent  to  (6)  is 

u(x,)  *  yei1lCX  +  x  <0 

(7) 

u(x)  =  eliuc  }  x  >  L 


Comparison  of  (6j  and  (7)  shows  u.iat  T  *  1/m  and  R=  u/m  •  The  solution  u(x) 
can  be  separated  into  ics  real  and  imaginary  parts  u^(x)  and  u. (x) .  Both 
ur(xJ  and  u^ixj  are  solutions  to  Schrodinger  equation  with 

urCxj  =  cos(kx;,  x>L 


uiU)  -  sin(kxj  ,  x>L. 

The  initial  conditions  for  the  solutions  u^(x)  an(j  u^(x)  at  x  =  L  are 
ur(Lj  =  cos(kU 
ur"iL)  *  -k  sin(kU 
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and 


=  sin(kU 

* 

U1(L)  =  k  cos (kL)  . 

Botn  solutions  can  be  projected  backwards  into  the  interval  (  -2’n’/^»  0)  • 

Using  the  orthogonality  of  the  sine  and  cosine  functions  in  this  interval,  the 
coefficients  U  ,  V  and  thus  R,  T  are  obtained. 

IV.  NUMERICAL  COMPARISON  OF  EXACT  AND  WKB  SOLUTIONS 

Consider  the  problem  of  a  plane  wave  incident  from  the  left  on  a  potential 
barrier  Vlx; .  As  shown  in  Figure  1,  V(x)  is  zero  for  x  <  0.  The  complex 
wave  amplitude  u(x)  satisfies  equation  (5).  In  terms  of  the  free  space  wave 
nranoer  k,  define 


pix)  =  Lk^  -V(x.)  ]  ,  for  x  <a 

qlxj  =  [V(x>  -k2jI/2}  for  x  >a 


wnere 


Via)  =  k2. 
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oo 


in 


ref  1 


Fig.  1  —  One-sided  scattering  potential  where  incident  and  reflected  waves  have 
free  space  wavenumber  k.  The  classical  turning  point  is  a. 


The  variable  a  is  caiiea  the  classical  turning  point.  According  to  the 
connection  rormuias  of  the  WKb  metnod,  solutions  to  the  left  and  the  right  of 
x  =  a  are  given  by 


uiUi  = 


cos 


p ( x) dx 


x  <  a 


U2U)  = 


yj  q  ( x ) 


exp 


[-/' 


q  UJdx 


Consider  tne  region  to  tne  left  of  x  =  a  and  set 


a 

p ( x) dx . 


19) 


The  function  u^x)  can  be  expressed  as 


=  uinu)  +  urefl(x) 


where  for  x  <  0,  tne  incident  and  reflected  waves  are 


u.  (x)  =  -  exp  -A)]  exp  (ikx) 

in  vr  4 


ure f  1  (x)  -  -  exp  [i(A-  -J-)  ]  exp  (-ikx) 

■s/TT" 

Thus  the  phase  shift  of  tne  reflected  wave  relative  to  the  incident  wave  at 


x  =  0  is 


0W  =  2A 


TT 


2  ' 
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predicted  by  the  WKB 


Ine  final  result  is  that  the  total  chase  shirt  O 

w 

method  is  tne  geometrical  optics  term  2A  minus  tt/2. 


The  exact  value  <j>  ^  Qf  the  phase  shift  can  be  obtained  from  the  approach 
described  in  section  III.  The  WKB  phase  smft  error 


E  =  -  ®e. 


has  been  computed  for  several  potentials.  The  first  two 
are  cases  wnere  the  WKB  method  is  rairiy  accurate.  The 
relatively  broad  quadratic  potential  barrier  of  height  E 
width  W  at  naif  energy  value  is  1U  free  space  wavelength 


examples  presented 
first  example  is  a 

m  =  16  .  The  barrier 
that  is 


W  =  10  X 

0 


where 


4  TT 


x2 

O 


j  E 
2  m 


The  potential  is  shown  in  Figure  2  and  the  phase  shift  error  is  given  Figure 
3.  It  can  De  seen  that  the  error  is  only  a  few  degrees  over  a  significant 
range  of  wave  numbers.  Since  the  quantity  A  defined  by  equation  (9) 
approaches  0  as  the  wavenumber  goes  to  zero,  the  WKB  predicted  phase  shift  is 
90°  for  infinitely  small  wavenumbers.  Since  the  true  phase  shift  approaches 
160°  for  small  wavenumbers,  the  limiting  value  of  the  phase  shift  error  is 
90°  as  the  wavenumber  goes  to  0.  This  limit  has  been  observed  in  this  and 
other  examples. 
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V(x) 


Fig.  2  —  Broad  quadratic  potential  barrier  with  half  energy  width  equal  to  10 
free  space  wavelengths  at  corresponding  wavenumber 


phase  error  (in  degrees) 


Fig.  3  —  Phase  shift  error  for  the  potential  barrier  of  Fig,  2 


Tne  next  example  cons  laersa  is  the  linear  potential 

vu;  =  x. 

As  the  wavelength  become  small,  it  can  be  seen  from  Figure  4  that  the  phase 
smft  error  goes  to  0.  This  can  be  expected  since  for  short  wavelengths  in  a 
linear  potential  the  assumed  conditions  for  the  WKB  method  and  the  connection 
tormuias  are  satisfied. 

beverai  examples  where  the  WKB  method  has  a  much  larger  error  will  now 
be  aescriQea.  The  first  example  is  a  simple  modification  of  tne  quadratic 
barrier.  The  quadratic  barrier  of  Figure  2  is  reduced  in  width  so  that 

W  -  X 

o 

wa  ere 

111.  =  I  E 

,2  2  m 

o 

Then,  as  is  shown  in  Figure  5,  the  pnase  shift  error  becomes  significantly 
larger.  The  increased  error  can  be  expected  since  the  potential  is  a  more 
rapidly  varying  function  of  position  and  also  the  exponentially  decayed  form 
of  tne  wave  tunc t ion  for  x  >  a  as  assumed  m  equation  (8)  becomes  a  poor 
approximation  for  narrow  barriers. 

Tne  potential  well  shown  in  Figure  b  has  unique  pnase  shift  properties. 
The  pnase  shift  error  as  given  m  Figure  7  exhibits  the  typical  trend  of 
becoming  smaller  with  decreased  wavelength  with  two  exceptions.  The 
exceptions  are  narrow  regions  where  the  phase  shift  error  is  an  extremely 
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phase  error  (in  degrees) 
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Fig.  5  —  Phase  shift  error  for  a  quadratic  barrier  with  half  energy  width  equal  to  1  free 
space  wavelength  at  corresponding  wave  number 
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Phase  shift  error  for  the  potential  of  Fig.  6 


rapiaiy  varying  function  of  traquency  and  are  associated  with  the  existence 
of  metastabie  or  resonance  states  in  the  terminology  of  Reference  [2]  .  As 
aiscussea  in  Reference  i2j,  the  condition  for  the  existence  of  a  resonance  is 


CIO) 


tie  re  x 


L  and  X2  are  the  classical  turning  points  at  the  two  sides  of  the 
potential  well  for  wmch 


V<-V  -  VU2; 


ihe  coacii cion  cor  resonance  as  specified  by  equation  (10)  is  again  based  on 
tne  WKB  me  tnoa  *  It  is  worthwnile  to  compare  the  accuracy  of  the  WKB 
prediction  against  an  exact  solution.  The  exact  solution  method  employed  is 
cased  on  a  t mite  element  scheme  applied  to  tne  Eigenvalue  problem 


— Si  +  V(x)u  =  Eu 
.  2 


dx 


using  Hermite  cuDic  polynomials  as  basis  functions.  here  VCx)  is  the  poten¬ 
tial  well  snown  in  Figure  b. 

Figure  7  has  two  "points  of  discontinuity 11  where  the  phase  shift  error  is 
essentially  a  discontinuous  function  of  wave  number.  Since  the  WKB  predicted 
phase  shift  is  for  this  example  a  relatively  slowly  varying  function  of  wave 
numDer,  one  would  expect  that  tnese  points  are  wave  numbers  for  which  the  true 
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pnase  shift  nas  a  ois  con  tinuity .  That  this  is  the  case  can  be  seen. from 
Figure  8,  wnich  snows  tne  true  phase  shift  in  a  neighbornood  of  the  first 

point  of  discontinuity.  Using  a  simixar  expansion  at  the  nigher  wave  numoer , 
these  points  are  determined  to  four  significant  figures.  The  table  shown 
below  compares  tne  predicted  with  the  actual  location  of  the  points  of 
discontinuity . 


k. 

l 

k2 

Actual 

i .  boa 

1.936 

WKb  Resonance  Condition 

1 .690 

1 .942 

Ei genv  aiue  Solution 

1.668 

1.933 

It  is  interesting  to  note  that  for  tne  higher  energies  and  wave  numbers  one 
can  expect  tne  WKB  condition  for  resonance  to  become  more  accurate  while 
tne  energy  eigenvalues  determined  a  numerical  scheme  tend  to  become  less 
accurate.  This  observation  is  consistent  with  results  presented  above. 

Tne  potential  obtained  from  the  first  hump  of  the  potential  well  is  shown 
in  Figure  y  ana  the  corresponding  phase  shift  error  curve  is  given  in  Figure 
10.  It  snould  be  noted  that  except  for  narrow  regions  aoout  tne  two  points  of 
discontinuity,  tne  curves  in  Figures  7  and  10  overlap.. 
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phase  error  (in  degrees) 


Fig.  8  —  Expansion  of  the  phase  shift  error  curve  of  Fig.  7  about  first  point  of  discontinuity 
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Phase  error  (in  degrees) 
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Fig.  10  —  Phase  shift  error  for  the  potential  barrier  of  Fig.  9 
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sided  problem  the  amplitude  of  the  reflection  coefficient  is  of  unit  magnitude 
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The  last  two  examples  presented  are  ledge  type  potentials.  The  ledge 
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Fig.  11  —  Ledge  type  potential  barrier  of  width  one  free  space  wavelength 
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Fig.  13  —  Phase  shift  error  for  a  ledge  type  potential  of  width  two  free  space  wavelengths 
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True  phase  shift  for  ledge  type  potentials 


and  the  reflection  process  is  cnaractenzed  by  the  phase  of  the  reflection 
coefficient.  The  numerical  approach  is  about  as  easy  to  implement  as  the  WKB 
approximation  and  eliminates  the  pnase  errors  typically  associated  with  the 
WKh  method.  One  novel  result  obtained  is  that  there  are  commonly  occurring 
index  of  refraction  functions  for  which  the  true  phase  shift  of  the  reflected 
wave  has  one  or  more  points  of  discontinuity.  The  implications  of  this  type 
of  pnase  response  on  communication  channels  can  be  explored  by  doing  numerical 
simulations  of  pulse  propagation  taxing  as  the  medium  transfer  function  one  of 
unit  amplitude  and  with  a  discontinuous  phase  term.  As  an  example,  the  lower 
point  ot  discontinuity  of  Figure  7  is  very  narrow  and  in  general  could  not  be 
expected  to  have  any  noticeable  effect  on  pulse  shape  while  the  upper  point  of 
discontinuity  is  sufficiently  broad  to  significantly  change  the  shape  of 
communication  signals.  For  HF  communication,  one  would  want  to  select  a 
signal  bandwidth  wnich  avoids  the  step  tunction  discontinuity  in  the  phase 
response  associated  with  ledges  in  the  election  density  profile.  As  can  be  / 
seen  from  Fig.  14,  this  effect  is  present  even  with  a  relatively  narrow  ledge 
corresponding  to  a  width  of  a  few  wavelengths. 

There  are  some  implications  of  tnis  study  to  the  plasma  inverse  problem. 
Since  the  phase  responses  of  the  potential  well  and  the  corresponding  single 
nump  potential  coincide  (see  Figures  7  and  10)  except  at  the  resonance  points, 
it  is  apparent  that  in  general  no  formal  solution  to  the  inverse  scattering 
problem  can  in  a  practical  sense  infer  the  existence  of  secondary  maxima  of 
electron  density  based  on  reflected  energy.  However,  in  the  special  case 
wnere 'resonance  states  occur,  the  existence  of  secondary  maxima  can  be 
investigated  by  sweeping  through  a  range  of  frequencies  with  a  narrow  band 
signal  and  recording  tne  phase  shift. 
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